Analysis of Collective Neutrino Flavor Transformation in Supernovae 
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We study the flavor evolution of a dense gas initially consisting of pure mono-energetic v e and Re- 
using adiabatic invariants and the special symmetry in such a system we are able to calculate the 
flavor evolution of the neutrino gas for the cases with slowly decreasing neutrino number densities. 
These calculations give new insights into the results of recent large-scale numerical simulations of 
neutrino flavor transformation in supernovae. For example, our calculations reveal the existence of 
what we term the "collective precession mode". Our analyses suggest that neutrinos which travel 
on intersecting trajectories subject to destructive quantum interference nevertheless can be in this 
mode. This mode can result in sharp transitions in the final energy-dependent neutrino survival 
probabilities across all trajectories, a feature seen in the numerical simulations. Moreover, this 
transition is qualitatively different for the normal and inverted neutrino mass hierarchies. Exploiting 
this difference, the neutrino signals from a future galactic supernova can potentially be used to 
determine the actual neutrino mass hierarchy. 

PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

In this paper we employ physical, analytic insights 
along with the results of large-scale numerical calcula- 
tions to study the nature of collective neutrino and an- 
tineutrino flavor transformation in supernovae. Although 
neutrino flavor transformation is an experimental fact, 
modeling this process in astrophysical settings can be 
problematic. In part, this is because nature produces 
environments where the number densities of neutrinos 
and/or antineutrinos can be very large. Examples of 
these include the early universe and environments as- 
sociated with compact-object mergers and gravitational 
collapse. In particular, core-collapse supernovae result 
in hot proto-neutron stars that emit neutrinos and an- 
tineutrinos copiously from the neutrino sphere. This im- 
plies inhomogeneous, anisotropic distributions for these 
particles. As their trajectories intersect above the proto- 
neutron star, their flavor evolution histories are quan- 
tum mechanically coupled [1]. The flavor content of the 
neutrino and antineutrino fields in and above a proto- 
neutron star will be a necessary ingredient for the in- 
terpretation of neutrino signals from a future supernova. 
It can also be an important, even crucial determinant of 
the composition of supernova ejecta [2] and possibly even 
the supernova explosion mechanism [3] . Consequently, if 
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we are to understand core-collapse supernovae, it follows 
that we must understand neutrino and antineutrino fla- 
vor evolution in them. 

Large neutrino number densities imply that neutrino- 
neutrino in addition to neutrino-electron forward scatter- 
ing sets the potential which governs neutrino flavor con- 
version. Because of the neutrino-neutrino forward scat- 
tering potential [4-6], neutrino flavor transformation in 
the early universe and near the supernova core can be 
very different from that in the vacuum or in an ordi- 
nary matter background only [7-17]. Recent analytical 
and numerical studies have revealed a new paradigm for 
neutrino flavor transformation in supernovae [18-26], one 
which is completely different from vacuum oscillations or 
the conventional Mikheyev-Smirnov-Wolfenstein (MSW) 
mechanism [27-29]. 

A particular aspect of this new paradigm is best dis- 
cussed in the following framework: The 2x2 neutrino 
flavor transformation problem can be described as the 
motion of isospins in flavor space, wherein Ve/v^^ and 
v e /Vfj,,T correspond to the "up" and "down" states of 
these isospins [21]. Both analytical and numerical stud- 
ies have suggested that a dense neutrino gas originally 
in a "bipolar configuration" (i.e., with the corresponding- 
neutrino flavor isospins or NFIS's forming two oppositely 
oriented groups) tends to stay in such a configuration 
even though each isospin group is composed of neutri- 
nos and/or antineutrinos with finite energy spread. In 
other words, neutrinos and antineutrinos with different 
energies can experience collective flavor transformation 
at high neutrino number densities. This is very different 
from the conventional MSW paradigm in which neutrinos 
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and antineutrinos with different energies undergo flavor 
transformation independently. 

A neutrino system with a bipolar configuration is also 
referred to as a "bipolar system" . Because supernova 
neutrinos are essentially in their flavor eigenstates when 
they leave the neutrino sphere, they naturally form a 
bipolar system. The neutrino sphere is in the very high 
density, electron degenerate environment near the neu- 
tron star surface. 

For a simple bipolar system consisting of mono- 
energetic v e and v e initially, it has been shown that the 
evolution of the system is equivalent to the motion of 
a (gyroscopic) pendulum [24]. Therefore, a bipolar sys- 
tem generally can evolve simultaneously in two different 
kinds of modes, i.e. the precession mode and the nu- 
tation mode, in analogy to the mechanical motion of a 
gyroscopic pendulum. In the extreme limit of large neu- 
trino number density, a bipolar system is reduced to a 
synchronized system, which is in a pure precession mode 
characterized by a common synchronization frequency 
[14]. The evolution of bipolar systems in the presence 
of an ordinary matter background has been studied in 
Rcfs. [21, 24] using corotating frames. 

Refs. [22, 23] have presented by far the most sophisti- 
cated, large-scale numerical simulations of neutrino fla- 
vor transformation in the coherent regime near the su- 
pernova core. For example, these simulations for the 
first time self-consistently treated the evolution of neutri- 
nos propagating along various intersecting trajectories. 
These simulations clearly show that the conventional 
MSW paradigm is invalid near the supernova core where 
neutrino fluxes are large. However, analytical models so 
far have only corroborated some of the features demon- 
strated by the simulations, and there are some obvious 
gaps between the analytical and numerical studies. 

One of the gaps is that current analytical models of 
bipolar systems assume constant neutrino number densi- 
ties, which is not true in supernovae. Using some simple 
numerical examples, Rcf. [24] has shown that some inter- 
esting phenomena observed in the simulations [22, 23] are 
related to varying neutrino number densities. For exam- 
ple, the energy averaged neutrino survival probabilities 
change as neutrino number densities decrease with the 
radius. In this paper we will show that if the neutrino 
number density decreases slowly as the system evolves 
out of the synchronized limit at high neutrino densities, 
the bipolar system will be dominantly in a precession 
mode. The neutrino flavor evolution seen in the numer- 
ical simulations is the combined effect of this precession 
mode and the nutation modes that are generated as a 
result of the finite rate of change in neutrino number 
densities. 

Another important gap between analytical and numer- 
ical studies is that most of the current analytical models 
assume homogeneity and isotropy of the neutrino gas, 
which is not true of the supernova environment. A re- 
cent analytical study which assumes an initial state of 
v e and v e with equal densities shows that the collec- 



tivity (referred to as "coherence" or "kinematic coher- 
ence" in Rcf. [24]) of the nutation modes tends to break 
down quickly among different neutrino trajectories in an 
inhomogeneous and anisotropic environment [26]. How- 
ever, the numerical simulations presented in Refs. [22, 23] 
employed more realistic supernova conditions where the 
initial v e and v e as well as v^ T and D^, T do not have 
the same number densities. These simulations do show 
some clear signs of collective flavor transformation. One 
important example is the hallmark pattern in the fi- 
nal energy-dependent neutrino survival probability which 
has a sharp transition energy across all neutrino trajec- 
tories (Fig. 3 of Rcf. [23]). We have further analyzed the 
numerical results obtained in the large-scale simulations 
mentioned above and found that, apart from the non- 
collcctive nutation modes, neutrinos propagating along 
different trajectories were in a single, collective preces- 
sion mode, ft is this precession mode that facilitates the 
mechanism suggested in Ref. [22] for producing the fore- 
mentioned hallmark pattern in the final neutrino survival 
probability. 

This paper is organized as follows. In Sec. II we will 
study the properties of a symmetric bipolar system ini- 
tially consisting of an equal number of v e and D e . We will 
use an adiabatic invariant of the system to obtain some 
analytical understanding of the evolution of such a sys- 
tem as neutrino number densities change. In Sec. Ill we 
will compare a simple asymmetric bipolar system with a 
gyroscopic pendulum. We will revisit the criteria deter- 
mining whether a bipolar system is in the synchronized 
or bipolar regime and clarify the description of bipolar 
oscillations. In Sec. IV we will show that an asymmet- 
ric bipolar system can stay roughly in a pure precession 
mode if neutrino number densities decrease slowly. We 
will also demonstrate some interesting properties of such 
a precession mode which can explain the results from the 
simple numerical examples of Ref. [24]. In Sec. V we will 
apply our simple analytical models to understand the nu- 
merical simulations presented in Refs. [22, 23] and offer 
some new analyses of these simulations. In Sec. VI we 
give our conclusions. 



II. SYMMETRIC BIPOLAR SYSTEM 

A. Flavor pendulum 

We start with a simple bipolar system initially con- 
sisting of mono-energetic v e and D e with energy E v and 
an equal number density n v . Throughout this paper we 
will assume 2x2 flavor mixing through the active-active 
channel. According to Ref. [21], the flavor evolution of 
a neutrino or antincutrino is equivalent to the motion 
of the corresponding neutrino flavor isospin, or NFIS, in 
flavor space. For a neutrino, the NFIS in the flavor basis 
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is defined as 



1 I 2Re(at e a„ T ) 
- 2Im(a* e cv) 

\KJ 2 - KJ 2 , 



(i) 



where a„ e and a Vr are the amplitudes for the neutrino 
to be in v e and another flavor state, say v T , respectively. 
For an antineutrino, the corresponding NFIS in the flavor 
basis is 



1 



s P = -- 



2Re(a j}<i a* /T ) 

2im(a Pe a*J 
Ian I 2 - |a p I 2 



(2) 



where a Pe and a Pr are the amplitudes for the antineutrino 
to be v e and v Tl respectively. 

To obtain a simple analytical understanding of collec- 
tive neutrino flavor transformation, we will assume, un- 
less otherwise stated, that the neutrino gas is homoge- 
neous and isotropic and that there is no ordinary mat- 
ter medium. Using the NFIS notation, the equations of 
motion (e.o.m.) for the NFIS's Si (neutrino) and s 2 (an- 
tineutrino) of this simple bipolar system are [21] 



— Si = Si X (^ V ,1 H V + «i> n iA2S 2 ), 

-^s 2 = s 2 x (fi v ,2H v + u„n„,isi), 



and the initial condition is 



where 



Si(0) = -S 2 (0) = ^, 



n V A = n„ 2 = n v 



(3a) 
(3b) 

(4) 
(5) 



are the number densities of neutrinos and antineutrinos, 
and e z is the unit vector in the z-direction in the flavor 
basis. Eq. (3) clearly shows that the motion of the NFIS's 
is similar to that of magnetic spins. In this "magnetic 
spin" picture, the "magnetic spins" Si and s 2 precess 
around a common "magnetic field" 



Hy = — e!j. sin 29 v + e f z cos 29 v 
with "gyro magnetic ratios" 

5m 2 



Mv,i 



"MV,2 = 



2E„ 



(7) 



At the same time, si and s 2 are also coupled to each 
other with a coefficient 



Mp 



-2y/2G F , 



(8) 



where Gf is the Fermi constant. 

In this paper we always take the squared difference of 
the two neutrino vacuum mass eigenvalues to be positive 
(8m 2 = m\ — m\ > 0) . Accordingly, the vacuum mixing 



angle 9 V varies within (0, 7r/2). A normal mass hierarchy 
corresponds to a mixing angle 9 V with < 9 V < 7r/4 and 
an inverted mass hierarchy corresponds to 7r/4 < 9 V < 
n/2. For an inverted mass hierarchy scenario, we follow 
Rcf. [24] to define 



7T 

2 



(9) 



which has < 9 V < 7r/4. We will loosely refer to both 9 V 
and 9 V as vacuum mixing angles. 1 

We first look at the scenario with n v being constant. 
With the definition of 



S± = n^isi ± n^, 2 s 2 , 
Eqs. (3) and (4) become [21] 



— S + = flyS- x H v , 

= MvS+ x H v + /i^S_ x S+, 

dr 

and 

S+(0) = 0, 



(0) = n v e\. 



(10) 

(11a) 
(lib) 



(12a) 
(12b) 



It is more convenient to work in the vacuum mass basis 
where the unit vectors are related to those in the flavor 
basis by 



e y= e y 



e x cos 29 v + e z sin 29 v , 
f 



Hv = — e„ sin 29 v 



e, cos26* v . 



(13a) 
(13b) 
(13c) 



Using Eqs. (11) and (12) one can check explicitly that 
vector S_ rotates in the eJJ-e^ plane while S + varies only 
along the e^-axis [21]. 

In reality neutrinos can experience collective oscilla- 
tions only if ti„ is large. The largeness of the neutrino 
number density in this simple bipolar system is naturally 



(6) measured by the ratio n„/n°, where 



5m 2 



4V2G F E„ ' 



(14) 



1 Refs. [21-23] have adopted a different convention for the inverted 
mass hierarchy scenario where 8 V defined here is the vacuum mix- 
ing angle and 6m 2 is taken to be negative. We note that these 
two conventions arc equivalent by the simultaneous transforma- 
tions \v\) <-> | V2) and \u T ) — > — \v T ). Correspondingly, one has 

Hv — > — Hv and , , — > — , , in flavor space. The x- and 

x \y ) x \y) 

y-components of a NFIS in the flavor basis under these two con- 
ventions are different by a sign for the inverted mass hierarchy 
case. 
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where 



FIG. 1: The equivalence of a pendulum and a symmetric 
{n v ,i = fii/,2 = n v ) bipolar system initially consisting of pure 
v e and v e . In the limit n„/n° S> 1 two NFIS blocks Si = 
n„,isi (for neutrino) and S2 = n„,2S2 (for antineutrino) are 
always roughly anti-aligned and Si points roughly in the same 
direction as the pendulum in flavor space. 



In the limit n„/n° 3> 1, the last term in Eq. (lib) dom- 
inates and 



d* 



S_ ~ (i v S_ x S+. 



(lib') 



Therefore, S_ roughly maintains a constant magnitude if 
the neutrino number density is large. As a result, si and 
S2 are always roughly anti-aligned although their direc- 
tions can be completely overturned in some scenarios. It 
is after this special property that "bipolar" flavor trans- 
formation was initially named [21]. 

We define d as the angle between vectors S_ and Hy, 
which varies within (0, ir) if S_ -e^ > and within (— ir, 0) 
if S_ • el < 0. We also define 



S + -el 



(15) 



With the initial condition in Eq. (12), we find that 
Eqs. (11a) and (lib') are equivalent to 



P-& 



sint?, 



P-d- 



(16a) 
(16b) 



Eq. (16) can be further reduced to a differential equa- 
tion of i3 of the second order: 



-ui 2 sini 



(17) 



2 Ref. [24] appears to have misunderstood the origin of the word 
"bipolar" by stating that the notation "bipolar oscillation" is a 



s/2G-p5m 2 n v 
El 



1/2 



(18a) 
(18b) 



is an intrinsic frequency of the system. Because Eq. (17) 
also describes the motion of a pendulum, we can view 
the flavor transformation in this simple bipolar system 
as the motion of a pendulum in the flavor space (Fig. 1). 
We note that the mass M of the "flavor pendulum" is 
irrelevant in this case. The only relevant parameter is 
the ratio between the magnitude of the acceleration field 
g and the length of the pendulum r, which is related to 
the intrinsic frequency by 

(19) 

The period of the flavor pendulum is (see, e.g., Ref. [30]) 
4 J Fs:(sin(^ max /2)) 



where 



T = 



K{k) EE f 

Jo 



r/2 



dC 



y/l - k 2 sin 2 C 



(20) 



(21) 



is the complete elliptic integral of the first kind [31]. 

The period of the simple symmetric bipolar system in 
Eq. (20) takes a simpler form if the vacuum mixing angle 
9 V or 6 V is small. For the normal mass hierarchy scenario 
with S v < 1, the pendulum motion is the same as that 
of a harmonic oscillator and 



rr, 2lT 



2tt 



^/2G-F5m 2 n lJ 



1/2 



(22) 



For the inverted mass hierarchy scenario with 9 V <C 1, 
expand Eq. (20) in terms of 8 V [31] and find that 



we 



41n(4/0 v ) 



UJ 



= 41n(4/0 v ) 



E v 



V / 2G , F<5m 2 n^ 



1/2 



(23) 



"misnomer 



The period of the bipolar oscillation in this limit has 
a logarithmic dependence on 8 V as pointed out in 
Ref. [24] . The periods of bipolar oscillations calculated 
using Eqs. (22) and (23) are in excellent agreement with 
the simple numerical examples in Ref. [21]. 

Ref. [24] has shown that the evolution of this simple 
bipolar system is equivalent to a pendulum motion for 
any n v (also see Sec. IIIB). In the limit n^/n" 3> 1 
the flavor pendulum described here is the same as that 
in Ref. [24] . This limit is of interest to analyses of the 
flavor evolution of supernova neutrinos and antineutrinos 
which have finite spread in their energy distributions, and 
therefore, may experience the collective flavor transfor- 
mation only when n v is large [21]. 
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B. Slowly varying neutrino number density 

If the neutrino number density n u varies with time, 
Eq. (3) is still valid but Eq. (11) is not. In this case, 
Eq. (16) is also valid as long as n„/n[J 3> 1. We note 
that $ and p$ comprise a canonically conjugate coordi- 
nate and momentum pair. In these variables the flavor 
pendulum has Hamiltonian 



H = nv 



cos d 



(24) 



Ref. [24] first noticed that the amplitude of flavor mix- 
ing in this bipolar system (or equivalently, the maxi- 
mal angular position i9 max of the flavor pendulum) de- 
creases with the neutrino number density n v . Drawing 
an analogy to the relation between the kinetic energy 
and the angular momentum of a pirouette performer, 
Ref. [24] suggested an intuitive explanation for this phe- 
nomenon: As n v becomes smaller, the effective mass 
m e a = /iy 1 ( 71 v / n v) m Eq. (24) increases. As a result, the 
kinetic energy of the flavor pendulum Ey m = p\ / (2m c g) 
is reduced with smaller n„ and the flavor pendulum can- 
not swing as high as before. 

We can arrive at the same conclusion for the scenarios 
n v/ n 1 3* 1 using Eq. (16). Let us compare the evolution 
of two flavor pendulums (a) and (b). We assume that the 
two pendulums have the same values of ("d,p$, n„) 
at instant t = to- We also assume that pendulum (a) has 
constant n v and that (b) has n v decreasing with time. 
After an infinitesimal interval At, both pendulums will 
have the same values of {•d,ps,p^) but pendulum (b) has 
smaller (d, n v ) than (a) does [see Eq. (16)]. This is equiv- 
alent to saying that both pendulums have the same an- 
gular position and potential well but pendulum (b) pos- 
sesses less kinetic energy £^in = #p#/2 than (a) does. As 
a result, pendulum (b) will not swing as high as (a) even 
if n v is constant for t > to + At. 

We note that neither Eq. (16) in this paper or Eq. (7) in 
Ref. [24] is equivalent to the original e.o.m. of the NFIS's 
[Eq. (3)] if n v is small and varies with time. Therefore, 
this explanation fails for n v /nt J; 1. 

To quantify the relation between the maximal angular 
position $ max of the flavor pendulum and the neutrino 
number density n v , we note that in the limit n„/n° 3> 1 



A 



Ps dti 



(25) 



is an adiabatic invariant of the pendulum motion (see, 
e.g., Ref. [30]). The integration in Eq. (25) is performed 
over one pendulum cycle with — i) max < i? <• ^max- 

The 

neutrino number density n v is taken to be constant dur- 
ing this cycle. 

If n v is constant, the Hamiltonian of the flavor pen- 
dulum is also a constant and is — /iv cosi9 max . Using 
Eq. (24) we obtain 



M 



'2n° 



- \J COS {> — COS t? n 



(26) 
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FIG. 2: (Color online) Flavor oscillations of the simple sym- 
metric bipolar system in a nearly adiabatic process [Eq. (34) 
with 7 = 10 and n v (0)/n° = 100] for 6 V = 0.01. The dashed 
line is the value of — cos #i as a function n„ computed from 
Eq. (3) with the initial conditions in Eq. (4). The solid line 
is — cos i? m ax computed from Eq. (30) . 



Combining Eqs. (25) and (26) we have 



2n° f° 

A=~i\ / \/cOS $ - C0S ''max dtf 

V n » -Aw 



= 16W^W(tf max )- 



(27a) 
(27b) 



Function W(£) in Eq. (27b) is defined as 

7^f o V cos (- cos I dC 
£(sin(£/2)) - cos 2 (^/2)K(sin(^/2)), (28b) 



w(0 



(28a) 



where 



,7T/2 

E(k) = / 

JO 



1 - k 2 sin 2 CdC 



(29) 



is the complete elliptic integral of the second kind [31]. If 
n^/n® 3> 1 and n„ varies slowly (adiabatic process), then 
i) max as a function of time satisfies the following relation: 



W(tf max ) W(20 v )i 



l n v (t) 
n„(0)' 



(30) 



An interesting scenario is that 6 V — tt/2 — 6 V <^ 1. In 
Hy and the probability for a neutrino 



this case, e z 
to be v e is 



P„ 



si ■ e z 



„f 1 — COS - !?! < 1 — COSfl 



(31) 
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where #1 is the angle between the directions of si and 
H v . Noting that W(26 v ) ~ W(n) = 1, we find that 



where 



(-Pi/ e i/ e )r 



0<P„ e „ e (t) < (P„ 



1 - COS W 



< n v {t) 
n„(0) 



(32) 



(33) 



is the maximal value that P„ c „ c may take for a given 
n v . In the above equation IU _1 (£) is the corresponding 
inverse function of W(£). 

For a concrete example, we assume that n v (t) has a 
linear dependence on time t: 



n v {t) = n„(0)(l - 7 Vv*), 



(34) 



where 7 is the adiabatic parameter. The larger the value 
of 7, the more adiabatic the process is. Taking 9 V = 0.01, 
we calculate the value of — cos$i as a function of n v for 
7 = 10 and n„(Q)/v° = 100 by solving Eq. (3) with the 
initial conditions in Eq. (4) . The results are plotted as the 
dashed line in Fig. 2. We have also computed $ max from 
Eq. (30) and plot — cos$ max as the solid line in Fig. 2. 
It is clear that — cos $ max outlines the upper envelope of 
— cos^i for n„/n° 3> 1. 

Using the analogy of harmonic oscillators, Ref. [24] has 
argued that, for the scenario with 9 V « 1, {Pv e v e )max 
should depend linearly on y/n^, at least when t9 max <C 1. 
However, it is clear that this conjecture is not true if i9 max 
is significant. In this case, (P„ e „ e ) max can be understood 
using the general form of W(£) and Eq. (33). On the 
other hand, we note that 



0„ 



(35) 



for $ max -C 1 [31]. Using Eqs. (30) and (35) we obtain 



■d 2 

max 



/ n v (t) 



(36) 



where to is an instant at which $ max (to) "C 1. Because 
P VeV , ~ (1 - cos<?i)/2 ~ tf?/4, we have 



{Pv e v e )max — 



<a*(*o) 



n v (t ) 



(37) 



for P VeVe <C 1. Therefore, (P Ve v e )max does depend lin- 
early on ^fnZ in the limit z? max <C 1. Note that this result 
only applies for n v /n® ^> 1. As mentioned above, our ar- 
gument about the adiabatic invariant fails for n„/n° < 1. 



III. ASYMMETRIC BIPOLAR SYSTEM 

A. Gyroscopic flavor pendulum 

We now consider a simple asymmetric bipolar sys- 
tem initially consisting of mono-energetic v e and v e with 



different but constant number densities. We note that 
Eq. (3) is still valid except that we now take n v \ = n„ 
and n v ,2 — an,, with a / 1 being a positive constant. 
Ref. [24] has shown that this asymmetric bipolar system 
is equivalent to a gyroscopic pendulum or a spinning top 
in flavor space for which 



S+(0) = 



1 - a 



-n v e{ 



S_(0) = ±±^n„e f ,. 



To see this we define 



Q^S_-^H V . 



(38a) 
(38b) 

(39) 



(Although we follow Ref. [24] in demonstrating the equiv- 
alence of an asymmetric bipolar system and a gyroscopic 
pendulum, we have adopted somewhat different nota- 
tions for our convenience.) Using Eqs. (lib) and (39) 
one sees that Q obeys the e.o.m. 



-^Q = M^Q x S+ 

and maintains a constant magnitude 



(40) 



Q = n v 



1 + a 



1 + 



+ 



1 + a 



cos26» v 



1 + a 



1/2 



(41) 



With the definition of 



_ s+ 



J = 

n v 

Q 

Eqs. (11a), (39) and (40) lead to 



r = ~Mv 



r x H v , 

n 



r x J. 



Using Eq. (43) one can easily show that 
a = J • r 



(42a) 
(42b) 

(43a) 
(43b) 

(44) 



is a constant of motion. From Eq. (43b) one obtains 

rxr = - Av (|)(,r-J). (45) 
We note that Eqs. (43a) and (45) arc equivalent to 



J = r x Mg, 

J = Mr x r + or, 



(46a) 
(46b) 
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where 



M = 



g = 4q h v- 



(47) 
(48) 



Therefore, this asymmetric bipolar system is indeed 
equivalent to a gyroscopic flavor pendulum. Specifically, 
r is the position vector of the bob, J is the total angular 
momentum, a is the internal angular momentum of the 
bob, M is the mass of the bob, and g is the acceleration 
field. The only difference between this pendulum and 
that shown in Fig. 1(a) is the spin of the bob. Hereafter 
we will loosely refer to both the symmetric and asym- 
metric bipolar systems as flavor pendulums. 

The motion of the gyroscopic flavor pendulum is the 
combination of a precession around Hy and a nutation 
with i9 m i n < z9 < i5 mm . Here is the polar angle of r with 
respect to = Hy, and z9 m i n and i9 max are the minimal 
and maximal values of $ during nutation. For the simple 
asymmetric bipolar system that we have discussed, one 
has $ max = | t=o • The value of i? m ; n can be determined 
as follows. Following Ref. [24] we define the total energy 
of the pendulum as 



Etot — E pot + -Ekin 

= Mg(l - H v • r) + 



2M 

, Mr 2 a 2 
= M, ? (l-c stf) + — + — . 



(49a) 
(49b) 

(49c) 



We note that E to t differs from the conserved total effec- 
tive energy of the NFlS's [21] by only a constant multi- 
plicative factor and a additive constant, and therefore, is 
also conserved. Because the motion of the pendulum is 
a pure precession around Hy when $ = i9 m i n , one has 

M5(costf min - cos?9 max ) = sin 2 tf m in^ 2 |tf=tf min , 

where ip is the azimuthal angle of r with respect to Hy. 
Using the conservation of the total angular momentum 
in the direction of the acceleration field g, one obtains 

er(cos tf min - cos # max ) = -M sin 2 tf min <^|^=^ min . (51) 

Combining Eqs. (50) and (51) we have 



COS iVin = ~T]+ + cos l?max) 2 + (1 ~ COS 2 tf ma x), 

(52) 

where 

^4W < 53 > 

B. Precession/nutation modes and 
synchronized/bipolar regimes 



sponding analogous flavor pendulum is undergoing pre- 
cession (nutation). The symmetric bipolar system dis- 
cussed in Sec. II corresponds to the limit a = 1 and is 
always in a pure nutation mode. In this limit a = 
and t? m in = 0, so the pendulum does not spin at all and 
simply swings in a fixed plane. Taking p$ = J • and 
d as the angle between r and Hy, one can obtain from 
Eq. (43) that [24] 



P-d = -Mv 



Q 



sin$ 



(54a) 
(54b) 



Eq. (54) is the exact version of Eq. (16). In the limit 
n^/n® » 1, Q ~ S_ and S_ approximately follows a 
plane pendulum motion as we have discussed in Sec. II A. 
If the neutrino number density n v is constant and the vac- 
uum mixing angle 6 V or 8 V is small, the bipolar systems in 
the pure nutation mode can experience almost complete 
flavor conversion during a nutation period. This is true 
for various initial configurations (see Table I in Ref. [21]). 

A bipolar system generally evolves simultaneously in 
both precession and nutation modes. However, if the neu- 
trino number density is large enough, it has been shown 
that a neutrino gas is in the synchronized mode with a 
characteristic frequency O sync independent of its initial 
configuration [14]. The criterion for synchronization can 
be written as [21] 



where 



|^S| > |(mv)|, 



S — ^ ' Tlv,iSi 



(55) 



(56) 



is the total NFIS, and 



(MV> = }^ — g2 = Q 



sync 



(57) 



is the average vacuum oscillation frequency The index i 
in Eqs. (56) and (57) denotes neutrinos or antincutrinos 
with a specific momentum. 

For the initial condition in Eq. (38) we note that, if 
n v is large, the total angular momentum of the flavor 
pendulum is dominated by its spin 



or 



(58) 



and 



We shall refer to neutrino flavor transformation as be- 
ing in the precession (nutation) mode when the corre- 



l-a 



a ~ 



(59) 
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In the limit 3 



x r, 1 + a n 
n„ > n v = 8— yzK, 



the parameter rj [see Eq. (53)] satisfies 77 ^> 1 and 



(60) 



COStL 



-v + v 1 



cosz9 n 



cos^ max . (61) 



As a result, the flavor pendulum roughly maintains a con- 
stant latitude and is essentially in the precession mode. 
One can explicitly show that in this case the flavor pen- 
dulum precesses around — Hy with a constant angular 
frequency [24] 



n ~ n 



1 



sync 



1 



(62) 



Therefore, a bipolar system is synchronized if n v 3> n* 
and the synchronized mode corresponds to a pure preces- 
sion mode with the synchronization frequency O sync as its 
precession frequency. We refer to the limit in Eq. (60) as 
the "synchronized regime" . We say that a bipolar system 
is in the "bipolar regime" if Eq. (60) is not satisfied. In 
this case it can be in both the precession and nutation 
modes. 

For the simple asymmetric bipolar system discussed 
here, Eq. (55) lead to 



n v > 2 



1 



(I -a) 



K 
4 



(63) 



Eqs. (60) and (63) differ by a constant multiplicative fac- 
tor. This reflects the fact that there is no sharp boundary 
between the synchronized and bipolar regimes. The sim- 
ple prescription for the synchronization frequency f2 sy nc 
in Eq. (57) allows a ready and practical application of 
the synchronization condition in Eq. (55) for neutrino 
and/or antineutrino gases with finite spreads in their en- 
ergy spectra. 

The ways in which the word "bipolar" has been used in 
the literature [21-24, 26] can be very confusing. There is 
a tendency to mistakenly identify the synchronized (bipo- 
lar) regime with the precession (nutation) mode. This is 
probably because a flavor pendulum can only precess in 
the synchronized regime and a symmetric bipolar sys- 
tem was once viewed as a typical bipolar system which 
is always in a nutation mode. However, the criterion 
determining whether a bipolar system is mostly in the 
precession or nutation mode is not the same as that for 
determining whether it is in the synchronized or bipolar 



3 Ref. [24] first noticed that the flavor pendulum with 8 V <C 1 
possesses little nutation for n„ 3> and is in the synchronized 
regime. We note that a flavor pendulum will have little nutation 
so long as Eq. (60) is satisfied. This result is independent of the 
value of 8 V . 



regime. A good example is that an asymmetric bipolar 
system can simultaneously be in both the precession and 
nutation modes in the bipolar regime. We also note that, 
while the precession frequency of the precession mode 
in the synchronized regime is determined from Eq. (57) 
and is independent of the neutrino number density n v , 
the precession frequency of a precession mode in the bipo- 
lar regime depends on n v (see Sec. IV C). 

The evolution of a bipolar system in the bipolar regime 
is also referred to as bimodal oscillations. Ref. [12] has 
shown that, in an asymmetric bipolar system initially 
consisting of mono-energetic v e and 9 e , the x- and y- 
componcnts of the polarization vectors of the neutrino 
and antineutrino (2si and — 2s2 in the NFIS notation) 
are bimodal as they are functions of two intrinsic peri- 
ods. It is clear that these two periods are related to the 
precession and nutation of the flavor pendulum. If v e and 
v e have different energies or the system starts with differ- 
ent neutrino/antineutrino species, one can demonstrate 
that such a system is equivalent to a flavor pendulum 
in some properly chosen corotating frame [21]. In this 
case the precession frequency f2 is shifted by the rotation 
frequency of the corotating frame. 



IV. PURE PRECESSION MODE OF 
ASYMMETRIC BIPOLAR SYSTEMS 

A. Pure precession mode 

Although a bipolar system tends to develop some nu- 
tation in addition to the precession mode in the bipolar 
regime, the actual mix of these modes depends on the ini- 
tial conditions as well as the system configuration. For 
example, a flavor pendulum precesses around — Hy with 
constant angular frequency SI without any nutation if 

Mgsmd = n,J 1 _ (64a) 
= n(Mfisini?cosi? + <7sini?) (64b) 

is satisfied, and the corresponding bipolar system is in 
the pure precession mode. For a gyroscope this is known 
as the "regular precession" . 

With varying neutrino number densities the problem 
is generally complicated. This is because almost all the 
parameters of the flavor pendulum (Q, M, g, a, etc.) de- 
pend on n v and the c.o.m. of a pendulum, Eq. (46), is not 
equivalent to that of the NFIS's if n v (t) is not constant. 
In this case, one has to use Eq. (3) to follow the evolu- 
tion of the bipolar system. Simple numerical examples 
presented in Ref. [24] seem to suggest that the evolution 
of a bipolar system with 6 V ~ it/ 2 can be dominantly in 
a precession mode after the system transitions from the 
synchronized regime into the bipolar regime. Here wc try 
to gain some analytical understanding of this precession 
mode by using the same simple bipolar system studied 
in Sec. Ill but with time- varying n v . 

We note that, in the synchronized regime (i.e., the 
limit of large n v ), both Si and S2 precess uniformly 
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FIG. 3: (Color online) Values of cos i?i (solid lines) and cos #2 (dashed lines) , where i9i and $2 are the polar angles of the NFIS's 
si and S2 with respect to Hv, as functions of neutrino number density n v for simple asymmetric bipolar systems. Panels (a), 
(b) and (d) have mixing angle 6 V — 0.01 and panel (c) has 9 V = 0.6. The asymmetry parameter is a = 0.8 for (a) and (c), 
and a = 0.2 for (b) and (d). The thin lines are computed from the original e.o.m. of NFIS's si and S2 [Eq. (3) with the initial 
conditions in (4)] assuming a nearly adiabatic process [Eq. (34) with n„(0)/n* = 2]. The adiabatic parameter is 7 = 40 for 
panels (a-c) and is 200 for panel (d). The thick lines are computed from Eq. (66) assuming that bipolar systems stay in the 
pure precession mode. The vertical dot-dashed lines in panels (a), (b) and (d) correspond to n v — n£. 



around Hy, and the motion of the NFIS's has a cylindri- 
cal symmetry around the axis along Hy. This symmetry 
is inherited from the e.o.m. of the NFIS's [Eq. (3)]. We 
consider an infinitely long process during which n„ is de- 
creased without preference to any azimuthal angle with 
respect to Hy- The cylindrical symmetry in the motion 
of the NFIS's around Hy is expected to be preserved in 
such a process, and si and S2 keep on precessing uni- 
formly around Hy without any wobbling. 

If this is true, vectors s 1; S2 and Hy must always be in 
the same plane, and Si and s 2 rotates around -Hy with 
the same angular frequency 



n 



1 



2sint9i 
1 



2sini? 2 



sin(i?i + 1? 2 ) 



(65a) 
(65b) 



where $1(2) is the angle between s 1 ( 2 ) and Hy. On 
the other hand, from Eq. (3) it can be shown that 
Hy • (si + CVS2) is time invariant even if n v changes with 
time. Consequently, we obtain the following two equa- 
tions for $1 and $ 2 : 

4sini9i sin ^2 = — a sini? 2 ) 

x sin(#i +# 2 ), (66a) 
(1 — a) cos2# v = cosz?i + a cos $2- (66b) 

We have solved Eq. (66) numerically for simple asym- 
metric bipolar systems with different choices of 9 V and 
asymmetry parameter a. The results are plotted in 
Fig. 3. For comparison, we have also solved numerically 
the original e.o.m. of the NFIS's, Eq. (3), for the same 
bipolar systems assuming that n v changes in the way 
described by Eq. (34). These results are also shown in 
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Fig. 3. Clearly, the polar angles $1 and $2 of the NFIS's 
S! and s 2 oscillate around those values determined from 
Eq. (66) as n v decreases. This is true not only for the 
bipolar systems with 6 V ~ tt/2 but also for those with 
other vacuum mixing angles. 

The results shown in Fig. 3 can be understood as fol- 
lows. Although Eqs. (3) and (46) are not equivalent over 
a long period for a time- varying n v , we may still view a 
bipolar system as a flavor pendulum over a short time 
interval during which n„ does not change much. Sup- 
pose that at instant t x the flavor pendulum precesses 
uniformly around Hy at latitude "&o(t{). In the adiabatic 
limit this precession continues as n v slowly changes, but 
the value of $0 changes with n v [i?o is a function of $1 and 
#2 which vary with n v according to Eq. (66)]. Of course, 
in realistic conditions, n v can only decrease with a finite 
rate, and the actual polar angle •& of the pendulum always 
"wobbles" (as a result of excitation of nutation modes) 
around $ with some nutation period T nut . However, if 
n v changes so slowly that 



z9 can be expected to closely follow tfo, and Eq. (66) be- 
comes an excellent approximation. This expectation can 
be verified by comparing panels (b) and (d) of Fig. 3 
where the evolution of two otherwise identical bipolar 
systems is calculated using different adiabatic parame- 
ters [7 = 40 in (b) and 200 in (d)]. With a much slower 
change in n„ for panel (d), the result obtained by solv- 
ing Eq. (66) becomes very close to that derived from the 
exact numerical calculations. 

The final values of $i| n „=o and ^2|n„=o of the bipolar 
system in the pure precession mode can be obtained as 
follows. The precession frequency £1 of the flavor pendu- 
lum cannot be 0, and therefore, cannot change sign as n v 
decreases. Because f2 = fl syn c > in the limit of large 
n v for a < 1 [Eq. (62)], we have Q > 0. Eq. (65b) would 
give > for n v = only if 

02k=O=7T. (68) 

This and Eq. (65) then give 

^|n„=o = Mv- (69) 

Combining Eqs. (68) and (66b) we obtain 

cosi9i|„ iy= o = (1 — a) cos2# v + a. (70) 

This agrees with the numerical results shown in Fig. 3. 
For the inverted mass hierarchy scenario with 6 V ~ 7r/2, 
we have e l ~ — Hy and 



VeV c \ n u =0 — 



1 — COS$i 



~ 1 — a, 



n„=0 



1 + COS $2 



f e y e \n v — — 



~ 0. 



(71a) 
(71b) 



However, we note that these results for n v = do not 
apply to realistic bipolar systems with finite spreads in 
the neutrino and antineutrino energy spectra as collective 
oscillations of these systems always break down before n„ 
reaches 0. 



B. Critical neutrino number density for the 
inverted mass hierarchy scenario with 8 V ~ tt/2 

In Fig. 3 one can see that, for the inverted mass hi- 
erarchy scenario with 9 V ~ tt/2, si and S2 begin to mis- 
align with Hy when n v is smaller than some critical value 
n£, and there seems to be discontinuity in d'i?i( 2 )/dn I/ at 
n v = n c v . To understand these results, we consider the 
limit where 8 V — tt/2. We define 



x 1(2) = sini? 1(2 ). 
For $1 ~ 7r and t? 2 — we have 



2 ' 



COS-(9i ~ — 1 -+ 

x 2 

COS t? 2 — 1 — -y • 



Combining Eqs. (66) and (73) we obtain 



(72) 

(73a) 
(73b) 



4XXX2 ~ (Xl — (XX2)— 7T 



X2 1 



x\ 



xx\\- 



Eq. (74) has the solution 
xl ~ ax\ ~ 2y/a 

where 



1 \in v <nl, 



4n° 



(74a) 
(74b) 

(75) 
(76) 



Therefore, for the limiting case with V = tt/2, S\ and S2 
only start to misalign with Hy when n v becomes smaller 
than rf v . 

One can also obtain the same value for n£ from the 
gyroscopic flavor pendulum analogy using Eq. (52). For 
$max = 7T) one always has $ m i n = it if 4 



4M 2 g 



> 1. 



(77) 



n„=0 



4 Ref. [24] pointed out the existence of the critical neutrino number 
density using this argument but gives n% ~ raj. 



11 




FIG. 4: (Color online) Scaled precession frequency 0//xv of 
the flavor pendulum as a function of neutrino number density 
n v . The pendulum is assumed to be in the pure precession 
mode. The vacuum mixing angle is 6 V = 0.01 for the solid and 
dashed lines, and 9 V = 0.6 for the dotted line. The asymmetry 
parameter is a = 0.8 for the solid and dotted lines and is 0.2 
for the dashed line. The vertical dot-dashed line and dot- 
dot-dashed line correspond to n„ = n£ for a — 0.8 and 0.2, 
respectively. 



Such a gyroscopic pendulum is known as a "sleeping top" 
because the pendulum "sleeps" in the upright position 
defying the effect of the gravity (see, e.g., Ref. [32] for 
more discussions). 

We note that the period of nutation T nut of the flavor 
pendulum is infinite if 9 V = ir/2. For a symmetric bipolar 
system T nut oc | \n9 v \ if 9 V ~ n/2 (see Sec. II A and also 
Ref. [24]). One expects similarly long nutation periods 
for asymmetric bipolar systems in the region where n v ~ 
n%. In the same region, d\ and $i change very quickly 
in the pure precession limit. As a result, the condition in 
Eq. (67) is usually violated in realistic environments and 
significant nutation can appear for n v < n c v [see, e.g., 
Fig. 3(b)]. 



C. Precession frequency 

Using Eq. (65) we have calculated the precession fre- 
quency f2 of the flavor pendulum for several scenarios 
assuming the pendulum is always in the pure preces- 
sion mode. The results are plotted in Fig. 4. The pre- 
cession frequency D, asymptotically approaches the syn- 
chronization frequency f^sync in the synchronized regime 
{riv/n* ^> 1) as n v becomes larger and larger. On the 
other hand, f2 changes steeply with n v in the bipolar 
regime (n v /n* < 1). As n„ reaches 0, il = fiy and the 
flavor pendulum precesses with the vacuum oscillation 
frequency of the dominant neutrino species {y e in this 
case). 

We also note that, for the inverted mass hierarchy sce- 



FIG. 5: (Color online) Ratio of the potential energy E pot of 
the flavor pendulum to its kinetic energy [Eq. (49)] as 

functions of neutrino number density n„ . The parameters and 
the meaning of the lines are the same as those in Fig. 4. 



nario with 9 V ~ 7r/2, the precession frequency Q reaches 
its maximum at n„ ~ rf v . Using Eqs. (65) and (75), we 
obtain 

fik=„c = 1 + -p^= for0 v = J, (78) 
1 - Vat 2 

which is 17.9 for a = 0.8 and 2.6 for a = 0.2. These 
values agree well with the solid and dashed lines in Fig. 4 
which assume 9 V — 0.01. 



D. Equipartition of energies? 

Ref. [24] observed that, for 9 V ~ 7r/2, the total energy 
of a flavor pendulum begins to be approximately equipar- 
titioned between its potential and kinetic energies when 
n v reaches n c v . Using Eq. (75) one can show that (see 
Appendix), for the extreme case 8 V = n/2, the ratio of 
energies is -E po t /-Ekin ~ 1 to 0(n c v jn v — 1) at n v ~ n c u if 
the flavor pendulum stays in the pure precession mode. 

Ref. [24] also mentioned "an important detail" that 
energy equipartition cannot hold all the way down to very 
small n v because -Epot has a finite positive minimum and 
Skin ultimately reaches 0. This is illustrated in Fig. 8 
of the same reference. However, we can show that (see 
Appendix) 



E 



pot. 



1 



(1 - a) cos 2 29 v 
1 + a — 2a cos 29 v 



if the flavor pendulum stays in the pure precession mode. 

In Fig. 5 we plot the ratio E vot / i^kin as a function n v 
for three different bipolar systems in the pure precession 
mode with various choices of a and 9 V . Indeed, for the 
bipolar systems with 9 V ~ 7t/2, -Epot/^kin reaches 1 at 
n v = n c v and does not change much for a significant range 
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of n v . This is especially true for a ~ 1. In the limit n„ = 
0, E pot /E kin ~ 0.94 and 0.5 for V ~ tt/2 but a = 0.8 and 
0.2, respectively. In the same limit, -Epot/^kin — 0.98 for 
9 V = 0.6 and a = 0.8. These results agree well with 
Eq. (79). 



V. NEUTRINO OSCILLATIONS IN 
SUPERNOVAE 

Refs. [22, 23] have presented two sets of simulations 
using the "single-angle" and "multi-angle" approxima- 
tions, respectively. These simulations together with the 
simple analytical and numerical models discussed in the 
previous sections represent approximations to the real 
supernova neutrino oscillation problem at three different 
levels of complexity. 

The analytical and numerical calculations performed in 
this paper assume that the neutrino gas is homogeneous 
and isotropic and is represented by two mono-energetic 
neutrino and/or antineutrino species. 

The single-angle simulations increase the complexity 
by allowing each neutrino species (4 in the 2x2 case) to 
have continuous energy distributions. It assumes that the 
flavor evolution histories of neutrinos propagating along 
different trajectories are the same as those of neutrinos 
emitted radially from the neutrino sphere. Although the 
single-angle approximation incorporates the angle depen- 
dence of neutrino-neutrino forward scattering into the 
"effective neutrino density" [22] , it still assumes that neu- 
trino on all trajectories evolve similarly. 

The multi-angle simulations are by far the most sophis- 
ticated treatment of the problem. In these calculations 
neutrinos and antineutrinos have not only continuous en- 
ergy distributions but also continuous angular distribu- 
tions. The most important improvement implemented in 
the multi-angle simulations is that the flavor evolution of 
neutrinos and antineutrinos (with a wide range of ener- 
gies) propagating along different trajectories is followed 
independently. 

In this section we will first apply our simple models to 
the single-angle calculations. We will discuss the onset of 
neutrino flavor conversion in both the inverted and nor- 
mal mass hierarchy scenarios (Sees. VA and VB). We 
will also investigate the precession mode of the neutrino 
gas in supernovae and its effects (Sec. VC). Finally, we 
will offer some new analyses of the multi-angle simula- 
tions and comment on the collectivity of neutrino flavor 
transformation in supernovae (Sec. VD). 



Onset of neutrino flavor conversion in the 
inverted mass hierarchy scenario 



start flavor conversion after the neutrino number den- 
sity n v drops below some critical value n£. Although 
the conclusion was made in the absence of an ordinary 
matter background, it has been shown that the evolu- 
tion of bipolar systems is not changed qualitatively even 
in a dominant matter background as long as 6 V ~ ir/2 
[21, 24]. 

For a rough estimate of the radius where n v = n£, we 
assume that the neutrino gas behaves in a way similar to 
the simple bipolar system initially consisting of v e and 
v e with energies E Ve = ff MeV and E Vc = 16McV, re- 
spectively. (These values are the same as the average 
energies of v e and v e in the simulations.) In a properly 
chosen corotating frame, the evolution of this simplified 
bipolar system is the same as that of a gas initial con- 
sisting of mono-energetic v e and v e with energy [21] 



E v 



~ 13 MeV. 



1 



(80a) 
(80b) 



With the luminosities of all neutrino species being the 
same and L v = 10 51 erg/s, the ratio of the number den- 
sities of the two neutrino species is 



L u /E Ue 
L v / Ev c 



0.69. 



(81) 



Therefore, the critical neutrino number density is [see 
Eq. (76)] 



-7^(1 -V5)" 

V2G F Ev 

6.24 x 10 28 cm" 3 



(82a) 
(82b) 



for a mass squared difference Sm 2 = 3 x lO- 3 eV 2 . Using 
the single-angle approximation we estimate the effective 
neutrino number density to be 5 



nf (r) = 



4nR*E v . 



~ (1.25 x lQ 61 cm~ 6 ) 




(83a) 



(83b) 



where R v = 11km is the radius of the neutrino sphere 
adopted in the simulation. Therefore, n° ff ~ n c v at radius 
r c ~ 52 km. 

From panels (c) and (d) of Fig. 8 in Ref. [22] one sees 
that, in the single angle simulation, the flavor conver- 
sion starts at the radius rx — 63 km. At r > rx the 



The simulations presented in Refs. [22, 23] for the in- 
verted neutrino mass hierarchy scenario all have V = 0.1 . 
According to the discussions in Sec. IV B, bipolar neu- 
trino systems with vacuum mixing angle 9 V ~ tt/2 can 



5 Eq. (83) is similar to Eq. (40) in Ref. [22] except that we here are 
not computing the net effective neutrino density and, therefore, 
ignore the contribution of antineutrinos. 
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z-components of NFIS's s„ e experience rapid oscillations 
which correspond to the nutation mode of the flavor pen- 
dulum. The estimated value of r c ~ 52 km and the ob- 
served value of rx — 63 km differ by ~ 10 km. This 
difference most likely arises because at r ~ r c the nuta- 
tion frequency T~ ut of the flavor pendulum is very small. 
Consequently, there is a delay before significant nuta- 
tion amplitude can develop. On the other hand, smaller 
nutation frequency implies less adiabatic evolution [see 
Eq. (67)]. So once developed, the nutation amplitude 
will be large. The oscillation amplitudes of (s VeZ ) and 
(sp e z) are indeed large as shown in Fig. 8 of Rcf. [22]. 

We note that the region (r c < r < rx) where the 
nutation modes are to be excited is roughly the same 
region where the chaos- like phenomenon shown in Fig. 12 
of Ref. [22] occurs. In this region the differences of two 
almost identical systems can grow exponentially as they 
evolve. 



B. Onset of neutrino flavor conversion in the 
normal mass hierarchy scenario 

The simulations presented in Refs. [22, 23] for the nor- 
mal neutrino mass hierarchy scenario all have 9 V = 0.1. 
According to the discussions in Sec. Ill, a bipolar neu- 
trino system with vacuum mixing angle fl, « 1 corre- 
sponds to a flavor pendulum that oscillates in a very lim- 
ited region near the bottom of the potential well, and 
therefore, does not experience much flavor transforma- 
tion. However, this conclusion only applies in the absence 
of an ordinary matter background. 

In the presence of a matter background, it has been 
shown that, in the synchronized regime, neutrinos and 
antincutrinos of all the species and energies go through 
an MSW-like resonance simultaneously in the same way 
as does a neutrino with the characteristic energy E sync 
in the conventional MSW picture [18]. A similar phe- 
nomenon may also occur in bipolar systems in the bipo- 
lar regime (i.e., outside the synchronized regime) as sug- 
gested in Ref. [21]. If this is true, the dominant neutrino 
species are changed from v e and v e to and D^. T . Us- 
ing the corotating frame, one can show that the evolu- 
tion of a Vp^-v^ T gas with the normal mass hierarchy is 
similar to that of a v e -v e gas with the inverted mass hi- 
erarchy [21], and the flavor pendulum is raised from the 
bottom position to the top position because of the change 
in the dominant neutrino species. Bipolar systems can 
subsequently develop nutation modes after the collective 
MSW-like resonance. 

From panels (a) and (b) of Fig. 8 in Ref. [22] one sees 
that, in the single angle simulation, the z-components of 
NFIS's s Ve suddenly change at radius rx — 88 km and os- 
cillate rapidly afterwards. This corresponds to the initial 
collective MSW-like resonance followed by the nutation 
modes. We note that the observed value of r x in the sim- 
ulation is larger than the value of 74 km estimated for the 
fully synchronized limit [22]. The difference arises partly 



because the MSW-like resonance actually occurs in the 
bipolar regime in this case. 



C. Precession mode and final neutrino survival 
probabilities 

As shown in Sec. Ill bipolar systems generally are in 
both precession and nutation modes. This is indeed seen 
in the single-angle simulations for both the normal and 
inverted mass hierarchies. In Fig. 8 of Ref. [22], the x- 
and y-components of NFIS's s„ e oscillate with an approx- 
imate phase difference of it/ 2, signifying precession in the 
x-y plane. 

In Sec. IV A we have argued that, in the absence of an 
ordinary matter background, the intrinsic precession an- 
gular velocity ft of the bipolar system as a whole should 
be in the same direction as that of a single neutrino of the 
dominant species. Therefore, we expect bipolar systems 
dominated by neutrinos to tend to precess around -Hy. 
In the presence of a matter background, bipolar systems 
will also tend to precess around the direction opposite to 

H e = -e { z V2G F n e , (84) 

where n e is the net electron number density. For the in- 
verted mass hierarchy with 6 V ~ ir/2 and H v — — e*, 
the intrinsic f2 of the flavor pendulum is roughly in the 
same direction as the precession stemming from the mat- 
ter background. In this case, the combined precession 
does not change direction. As a result, the precession 
of NFIS's is always roughly around — Hy — e 2 for the 
inverted mass hierarchy. 

For the normal mass hierarchy with 8 V <C 1 and 
Hy — e z , the intrinsic $7 of the flavor pendulum is 
roughly in the opposite direction to that of the precession 
due to the matter background, and the combined preces- 
sion may change its direction. However, it is expected 
that the matter background becomes negligible after the 
collective MSW-like resonance. Therefore, the NFIS's 
precess roughly in the direction of — Hy ~ — e 2 in the 
region r > rx for the normal mass hierarchy. According 
to Fig. 8 of Ref. [22] NFIS's indeed precess around e z 
for the inverted mass hierarchy scenario and —e z for the 
normal mass hierarchy scenario. 

Ref. [22] has shown that (see Fig. 9 of that reference), 
for the inverted mass hierarchy scenario and at large 
radius, neutrinos with energies below Eq — 9 MeV are 
mostly in their initial flavors while neutrinos with larger 
energies and most antineutrinos can be completely con- 
verted to other flavors in the limit of large neutrino lumi- 
nosity L v . For the normal mass hierarchy scenario, neu- 
trinos with energies below Eq — 9 MeV are almost com- 
pletely converted to other flavors while neutrinos with 
larger energies and nearly all antineutrinos are mostly 
in their initial flavors in the limit of large L v . Ref. [22] 
suggested that this phenomenon is related to the preces- 
sion of NFIS's when neutrino number densities decrease 
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and the bipolar configuration starts to break down (see 
Fig. 10 of that reference). 

We note that the precession of NFIS's due to the mat- 
ter background H e is the same for all neutrinos, and we 
can essentially ignore it in a reference frame T\ rotating 
with angular velocity H e . In this corotating frame T\, 
the e.o.m. of NFIS s, is 



dt 



Sj x H 



tot 



= S l X (jUy^Hy + H„) 

= S t X (yUy.iHv + flv ^ n v,3 s j) 



(85a) 
(85b) 
(85c) 



where H„ is the effective "magnetic field" generated by 
all other NFIS's. We assume that all NFIS's and H„ ro- 
tate with a constant angular velocity -fiHy. The prob- 
lem becomes very simple in a reference frame Ti which 
rotates relative to T\ with angular velocity — £!H V - In 
T2 



dt 



Si = Si x H 



tot 



= s,x [(/i V ,i - 0)H V + H„], 



(86a) 
(86b) 



where both H to t and H„ are not rotating 6 . 

We first look at a NFIS Sj corresponding to a neutrino 
which is initially pure v e at the neutrino sphere. Be- 
cause \x v = — 2v / 2Gf is negative and the neutrino gas is 
initially dominated by v e , the NFIS Sj must be roughly 
antialigned with H to t — H„ when neutrino number den- 
sities are large. As neutrino number densities decrease 
to 0, Htot — * (Mv.i — 0)Hy- If neutrino number densi- 
ties decrease so slowly that the process is adiabatic, Sj 
will stay antialigned with H tot . At n v = the NFIS 
Si can be either aligned or antialigned with Hy depend- 
ing on whether /ly.i is smaller or larger than Q. For 
the inverted mass hierarchy scenario with 6 V ~ tt/2 and 
Hy — — e*, Si is roughly aligned with e f z ~ — Hy if 
/iy = Sm 2 1 '(2E V ) > fl and antialigned with e z otherwise. 
Accordingly, at large radii neutrinos starting as v e are 
still mostly in the u e flavor if their energies are below 



Eq = 



8m 2 
~27T 



(87) 



and are almost completely converted to other flavors oth- 
erwise. In words, one has 



P V „{E V ) ~ 



1 \iE v < E c , 
if E v > E c . 



(88) 



We have ignored the rotation of Hy in the corotating frames T\ 
and Ti because Hy is roughly in the same or opposite direction 
as the rotation axis for 6 V -C 1 or 8 V ~ n/2. 



TABLE I: Final neutrino survival probabilities in the adia- 
batic limit for large L„ in both the normal and inverted neu- 
trino mass hierarchy scenarios. 





9 V <S 1 (normal) 


9 V ~ 7r/2 (inverted) 


P VV {E V < E c ) 





1 


P VV {E V > E c ) 


1 





Pvv{Ev) 


1 






One can estimate the final neutrino survival probabilities 
for other cases in a similar fashion. We have summarized 
the results for the relevant scenarios in Table I. 

In this analysis we have assumed to be constant. 
This analysis is expected to hold as long as the process 
is more or less adiabatic and neutrino number densities 
decrease slowly. The predictions from this simple anal- 
ysis generally agree with the results of single-angle nu- 
merical simulations presented in Fig. 9 of Ref. [22] . The 
agreement is especially good for large neutrino luminosi- 
ties and in the neutrino sector for which P UV {E U ) has a 
relatively sharp transition or jump at E v ~ Eq. This 
pattern can be taken as a hallmark of collective neutrino 
flavor transformation because it results from a neutrino 
background that is in a collective precession mode. 



D. Collectivity and non-collectivity of neutrino 
oscillations in supernovae 

Single-angle simulations assumed that neutrinos of the 
same species and energy all evolve in the same way even 
if they are emitted in different angles from the neutrino 
sphere. This is not necessarily always a good approxima- 
tion as neutrino-neutrino forward scattering is angle de- 
pendent and the neutrino density distributions in the su- 
pernova environment are inhomogeneous and anisotropic. 
Even if neutrinos moving along various trajectories all 
have the same flavor content at a given radius, they pro- 
duce different refractive indices for neutrinos propagating 
in different directions. Therefore, neutrinos propagat- 
ing in different directions are expected to have different 
flavor evolution histories. On the other hand, neutri- 
nos propagating along different trajectories are coupled 
to each other through neutrino-neutrino forward scat- 
tering. The correlations among different neutrino tra- 
jectories are especially strong when neutrino fluxes are 
large. The inhomogeneity/anisotropy of the environment 
and the correlation among different neutrino trajectories 
act as two opposite "forces" which try to break and up- 
hold, respectively, the collective aspect of neutrino os- 
cillations in supernovae. At the moment it is difficult 
to perform an analytical study that can clearly predict 
which force wins. Our limited goal here is to gain insight 
into the issue of collectivity of neutrino flavor transforma- 
tion by analyzing the multi-angle simulations presented 
in Refs. [22, 23]. 
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FIG. 6: (Color online) The evolution of ^-components of the average NFIS's of supernova neutrinos in the flavor basis obtained 
from the multi-angle simulations as presented in Fig. 4 of Ref. [22]. The left panels are for neutrinos emitted as v B , and the right 
panels are for antineutrinos emitted as v e . The top panels employ the normal mass hierarchy with 8 V — 0.1, and the bottom 
panels employ the inverted mass hierarchy with 9 V = 0.1. The solid, dotted and dashed lines give (sl, x ) along the trajectories 
with cosOo = 1, 0.5 and 0, respectively, where Go is the emission angle defined in Ref. [22]. The neutrino mass-squared 
difference is taken to be Sm 2 = 3x 10~ 3 eV 2 . 



As shown in Fig. 2 of Ref. [23] , the flavor evolution of 
neutrinos on each trajectory in multi-angle simulations 
looks qualitatively similar to that in the corresponding 
single-simulations. However, the oscillations in neutrino 
survival probabilities P vv and P pp have different frequen- 
cies for different trajectories. For vacuum angle 6 V -C 1 or 
6 V ~ 7r/2, the oscillations in P vv and Pi,v represent the 
nutation of the flavor pendulum. Therefore, the nuta- 
tion modes of neutrinos propagating along different tra- 
jectories cannot be viewed as collective. Indeed, it has 
recently been shown that the nutation modes for sym- 
metric bipolar systems can quickly develop large phase 
differences and "de-cohere" for neutrinos propagating in 
different directions [26]. 

On the other hand, Fig. 3 of Ref. [23] shows that 
Pw{Ev) obtained by multi- angle simulations has a sharp 
transition at E v = Ec as in single-angle simulations. In 
addition, the value of Eq is approximately independent 



of neutrino trajectory direction. If this transition is re- 
lated to the precession mode of neutrinos as suggested 
by Ref. [22] and further explained here, NFIS's corre- 
sponding to neutrinos propagating in different directions 
must precess with the same frequency. In Fig. 6 we have 
plotted {s f vx ), the ^-component of the average NFIS's in 
the flavor basis, as functions of radius r for three repre- 
sentative trajectories obtained from the multi-angle sim- 
ulations of Refs. [22, 23]. One indeed observes that the 
NFIS's along various trajectories are approximately in a 
single collective precession mode and precess around ±e z 
with approximately the same frequency. 



VI. CONCLUSIONS 

We have investigated the simple symmetric bipolar sys- 
tem using the flavor pendulum analogy. We have shown 
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that an adiabatic invariant of the pendulum motion can 
be used to study the evolution of such a bipolar system 
when neutrino number densities change slowly with time. 
We have also studied an asymmetric bipolar system us- 
ing the gyroscopic pendulum analogy. As a gyroscopic 
pendulum, a bipolar system generally can be in both the 
precession and nutation modes simultaneously except in 
the synchronized regime where only precession is possi- 
ble. 

We have shown that an asymmetric bipolar system can 
stay mostly in a pure precession mode as it transitions 
from the synchronized regime into the bipolar regime if 
neutrino number densities decrease slowly. The preces- 
sion frequency of the system generally varies with the 
neutrino number density and approaches the synchro- 
nization frequency in the synchronized regime. For the 
inverted mass hierarchy case with mixing angle 8 V ~ n/2, 
we have calculated a more accurate value of the criti- 
cal neutrino number density below which bipolar systems 
can start flavor transformation. Because supernova neu- 
trinos naturally form asymmetric bipolar systems, these 
analyses could be useful for understanding the qualitative 
features of neutrino oscillations in supcrnovae. 

We have further analyzed the recent numerical sim- 
ulations of neutrino oscillations in supernovae. These 
large-scale simulations suggest that neutrinos traveling 
on intersecting trajectories and experiencing destructive 
quantum interference nevertheless can be in the collec- 
tive precession mode. This mode can result in sharp 
transitions in the final energy-dependent neutrino sur- 
vival probabilities P VV (E V ) across all trajectories. This 
sharp transition in P vl/ (E v ) can be taken as a hallmark of 
collective neutrino flavor transformation. Moreover, this 
transition occurs differently for the normal and inverted 
neutrino mass hierarchies. Based on this difference, the 



neutrino signals from a future galactic supernova poten- 
tially can be used to determine the actual neutrino mass 
hierarchy. 



APPENDIX: POTENTIAL AND KINETIC 
ENERGIES OF AN ASYMMETRIC FLAVOR 
PENDULUM 

Let us compare the kinetic energy Skin of the flavor 
pendulum with its potential energy E po t [see Eq. (49)] 
in the pure precession mode. The potential energy is 
defined as 



E pot = Mg(l - H v ■ r) = ^ v (q - q • H v ), 
where 

Q ,n° 
q = — = Si - as 2 H H v . 

Its kinetic energy is defined as 



E 



kin 



where 



Mv f nv \ J 2 



2M 



J = Si + as 2 . 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



According to Eq. (75), we have sin 2 z9i ~ a sin 2 $2 and 
sin 2 d 2 ~ 26/ if 



S = -a- - 1 « 1. 

n v 

It is straightforward to show that 



(A.5) 



J 




1 + a 2 - 2acos(tfi + i? 2 ) (ri 



I^-^ + e^ 2 ). 



+ — (costfi - acosi? 2 ) + — 



(A.6a) 

(A.6b) 
(A.6c) 



In addition, we have 



Combining Eqs. (A. 6c) and (A. 7b), we obtain 



q • H v = si • H v - as 2 ■ H v + 

1 + a n° _ 

~ — + — + V^s. 

2 n v 



(A.7a) 
(A.7b) 



Spot — MV 

= Mv 



1 + a 



(1 + V^) 2 



1 + a 



(A.8a) 
(A.8b) 
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Similarly, we have 

j2 _ l + a 2 + 2Q:cos(#i+fl 2 ) 



and 



1 + 5 



(A.9a) 
(A.9b) 

(A.lOa) 
(A.lOb) 



Therefore, we obtain £ po t/Ski n — 1 to 0{8) in the limit 
<f«l. 



One can also estimate the potential and kinetic ener- 
gies of the flavor pendulum in the limit n„/n° < 1. In 
this limit we have si • Hy — cosi9i/2, s 2 • Hy — —1/2, 
and si • s 2 ~ — cosi?i/4. Therefore, 



tf = -ns,)2- :.'( (8i-Hv- 082-H V )+( J 



(6v 1+2 (l) (Sl ' Hv - aS2 ' Hv)+ (l) (Sl - aS2)2 

^ + si • H v - «s 2 • H v + I (^r) [(si - as 2 ) 2 - ( Sl • H v - as 2 • H v ) 2 



~ q Hy + 



n„ \ sin 2 i?i 



?J 8 ' 



(A.lla) 

(A.llb) 
(A.llc) 
(A.lld) 



which gives 
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Mv / n v , . 2 



(A.12) 



Similary, we have 

Ekin = f^(^y sl+aS2 f (A.13a) 

~ ^ (l + a 2 - 2acos<?i). (A.13b) 

Combining Eqs. (A.12) and (A. 13b), we obtain 



pot 



sin 2 i?i 



l + a 2 — 2a cos 



(A.14) 



In the limit n„/nJJ <C 1, cosi9i ~ (1 — a) cos20 v + a, and 
the above equation reduces to 



Sp 0t _ _ (1 - a) cos 2 26> v 
_E kin — 1 + a - 2a cos 29 v ' 



(A.15) 



This work was supported in part by NSF grant PHY- 
04-00359, the TSI collaboration's DOE SciDAC grant at 
UCSD, and DOE grant DE-FG02-87ER40328 at UMN. 
This work was also supported in part by the LDRD Pro- 
gram and Open Supercomputing at LANL, and by the 
National Energy Research Scientific Computing Center 
through the TSI collaboration using Bassi, and the San 
Diego Supercomputer Center through the Academic As- 
sociates Program using DataStar. We would like to thank 
A. Fricdland, T. Goldman, J. Hidaka and M. Patel for 
valuable conversations. 



[1] Y. Z. Qian and G. M. Fuller, Phys. Rev. D51, 1479 [2] Y.-Z. Qian, G. M. Fuller, G. J. Mathews, R. W. Mayle, 
(1995), astro-ph/9406073. J. R. Wilson, and S. E. Woosley, Phys. Rev. Lett. 71, 



18 



1965 (1993). 

[3] G. M. Fuller, R. W. Mayle, B. S. Meyer, and J. R. Wilson, 

Astrophys. J. 389, 517 (1992). 
[4] G. M. Fuller, R. W. Mayle, J. R. Wilson, and D. N. 

Schramm, Astrophys. J. 322, 795 (1987). 
[5] J. T. Pantaleone, Phys. Rev. D46, 510 (1992). 
[6] G. Sigl and G. Raffelt, Nucl. Phys. B406, 423 (1993). 
[7] S. Samuel, Phys. Rev. D48, 1462 (1993). 
[8] V. A. Kostelecky and S. Samuel, Phys. Rev. D49, 1740 

(1994) . 

[9] V. A. Kostelecky, J. T. Pantaleone, and S. Samuel, Phys. 
Lett. B315, 46 (1993). 
[10] V. A. Kostelecky and S. Samuel, Phys. Lett. B318, 127 
(1993). 

[11] V. A. Kostelecky and S. Samuel, Phys. Rev. D52, 621 

(1995) , hep-ph/9506262. 

[12] S. Samuel, Phys. Rev. D53, 5382 (1996), hep- 
ph/9604341. 

[13] V. A. Kostelecky and S. Samuel, Phys. Lett. B385, 159 

(1996) , hep-ph/9610399. 

[14] S. Pastor, G. G. Raffelt, and D. V. Semikoz, Phys. Rev. 

D65, 053011 (2002), hep-ph/0109035. 
[15] A. D. Dolgov et al., Nucl. Phys. B632, 363 (2002), hep- 

ph/0201287. 

[16] Y. Y. Y. Wong, Phys. Rev. D66, 025015 (2002), hep- 
ph/0203180. 

[17] K. N. Abazajian, J. F. Beacom, and N. F. Bell, Phys. 
Rev. D66, 013008 (2002), astro-ph/0203442. 



[18] S. Pastor and G. Raffelt, Phys. Rev. Lett. 89, 191101 

(2002), astro-ph/0207281. 
[19] A. B. Balantekin and H. Yiiksel, New J. Phys. 7, 51 

(2005) , astro-ph/0411159. 

[20] G. M. Fuller and Y.-Z. Qian, Phys. Rev. D73, 023004 

(2006) , astro-ph/0505240. 

[21] H. Duan, G. M. Fuller, and Y.-Z. Qian, Phys. Rev. D74, 

123004 (2005), astro-ph/0511275. 
[22] H. Duan, G. M. Fuller, J. Carlson, and Y.-Z. Qian, Phys. 

Rev. D74, 105014 (2006), astro-ph/0606616. 
[23] H. Duan, G. M. Fuller, J. Carlson, and Y.-Z. Qian, Phys. 

Rev. Lett. 97, 241101 (2006), astro-ph/0608050. 
[24] S. Hannestad, G. G. Raffelt, G. Sigl, and Y. Y. Y. Wong, 

Phys. Rev. D74, 105010 (2006), astro-ph/0608695. 
[25] A. B. Balantekin and Y. Pehlivan, J. Phys. G34, 47 

(2007) , astro- ph/0607527. 

[26] G. G. Raffelt and G. G. R. Sigl (2007), hep-ph/0701182. 

[27] L. Wolfenstein, Phys. Rev. D17, 2369 (1978). 

[28] L. Wolfenstein, Phys. Rev. D20, 2634 (1979). 

[29] S. P. Mikheyev and A. Y. Smirnov, Yad. Fiz. 42, 1441 
(1985), [Sov. J. Nucl. Phys. 42, 913 (1985)]. 

[30] L. D. Landau and E. M. Lifshitz, Mechanics (Pergamon, 
Oxford, 1976), 3rd ed. 

[31] I. S. Gradshteyn and I. M. Ryzhik, Table of integrals, 
series, and products (Academic, Boston, 1994), 5th ed. 

[32] J. B. Scarborough, The Gyroscope, Theory and Applica- 
tion (Interscience, New York, 1958). 



